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Abstract 

We investigate the internal structure of charged black holes with a spherically symmetric model in 
the double-null coordinate system. Hawking radiation is considered using the S-wave approximation 
of semiclassical back-reaction and discharge is simulated by supplying an oppositely charged matter 
to the black hole. In the stage of formation, the internal structure is determined by the mass and 
charge of collapsing matter. When the charge-mass ratio is small, a wormhole-like internal structure 
is observed. The structure becomes analogous to the static limit as the ratio reaches unity. After the 
formation, mass inflation induces large curvature in the internal structure, which makes the structure 
insensitive to the late-time perturbations. The internal structure determined from the formation 
seems to be maintained during a substantial mass reduction. The discharge and neutralization of 
charged black holes is also investigated for both non-evaporating and evaporating cases. Finally, we 
discuss the implications of the wormhole-like structure inside of charged black holes. 
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Figure 1: Penrose diagrams of static charged black holes for (a) M > Q, (b) M = Q, and (c) M < Q. 

1 Introduction 

Charged black holes have complicated internal structures. The Reissner-Nordstrom (RN) solution, a static 
electro-vacuum solution with spherical symmetry, is a good starting point from which to look over the 
structure of charged black holes pQ [2] . The RN solution has three different types of structure depending 
on the mass (M) and charge (Q) of the black hole (Figure [Q. All of them have time-like singularities, 
and they induce Cauchy horizons which are boundaries of the causally undetermined regions |3j. 

The existence of the time-like singularity and the Cauchy horizon has caused some issues related 
to cosmic censorship conjecture. If Q becomes greater than M, the naked singularity will violate weak 
cosmic censorship [3]. However, it is known that the electrostatic repulsion prevents such configurations 
[5]. On the other hand, if an observer can pass through the Cauchy horizon, he or she will see the effects 
of the singularity violating strong cosmic censorship However, in realistic charged black holes, the 
Cauchy horizon becomes a null curvature singularity via mass inflation [6] \J, [2] [9] . It plays a role as an 
impenetrable barrier to the observer, preserving strong cosmic censorship [7] [3]. The back reaction of mass 
inflation also deforms the causal structure of charged black holes, as shown in Figure [3J 

As time passes, M and Q evolve via Hawking radiation and discharge, which can deform the structure 
of black holes. When M approaches to Q, a non-extreme charged black hole will approach to an extreme 
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Figure 2: Charged black holes with mass inflation scenario. The collapsing matter generates mass in- 
flation which induces the space-like singularity, and the inner-Cauchy horizon becomes a null curvature 
singularity. 

one. On the other hand, when Q is decreased to nearly zero, the charged black hole will become a neutral 
one. In both cases, proper descriptions of the transition processes are required. 

To figure out the dynamic process of charged black holes with quantum effects, some numerical ap- 
proaches have been studied [5] [TU] [U [H] [E] Q3] [E] ■ When Hawking radiation was considered, the sepa- 
ration between the inner horizon and the Cauchy horizon was observed [121. The inner horizon bends in 
a space-like direction and approaches to the time-like outer horizon. This explains how a non-extreme 
charged black hole approaches to an extreme one where the two horizons coincide [T2 [15 . On the other 
hand, as a charged black hole takes in some oppositely charged matter, its space-like inner horizon shrinks 
and becomes a singularity [15] . This corresponds to the space-like singularity of a neutral black hole, and 
the Cauchy horizon from this singularity eliminates the internal structure. This explains how the structure 
of charged black holes becomes a neutral structure. According to the points mentioned above, we obtain 
the causal structure of charged black holes in Figure [3] [15 . 

In this paper, we focus on the internal structure of charged black holes. We examined the horizon 
structure inside of the outer horizon for important stages in time evolution: formation, evaporation, and 
neutralization of the black hole. Therefore, the present work is also a study of the time evolution of 
charged black holes. In Section [2j the numerical scheme is introduced. In Section [3j we briefly review the 
time evolution of charged black holes. In Section^ the formation stage of charged black holes with various 
M and Q is studied. In Section [5J we investigate the response of charged black holes due to a substantial 
mass reduction using an exotic matter field. In Section [6l the neutralization process is discussed in detail. 
In Section we explore the wormhole-like structure inside of charged black holes. The convergence and 
consistency of the results are discussed in Appendix [A"l 
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Figure 3: Causal structure of dynamical charged black holes. Evaporation and neutralization are consid- 
ered. 
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2 Model for a dynamical charged black hole 



2.1 Equations and integration schemes 

We followed the numerical setup of the previous researches [5] [10] [11] [12] [13] [15] [16] [17] . In this paper, we 
set G = c = fee = 4-7reo = 1, and mpi = Zpi = vE> where mpi and /pi are the Planck mass and length, 
respectively. 

We described a gauge-invariant Lagrangian for a complex massless scalar field <j) coupled with an 
electromagnetic field A a , and an exotic real scalar field ip [2] 

C = - (<f> ]a + ieA a <P)g ab (4,. b - ieA b ~4>) + ^ a g a % b - ^F ab F a \ (1) 

where F ab = A b . M — A a - b and e is the gauge coupling. We also considered the semiclassical Einstein equation 
to include gravity and Hawking radiation, 

G ab = Stt (Tg + (T«)) , (2) 

where (T5) is the renormalized energy-momentum tensor. 

Simply put, we make various charged black holes with (f> and observe the influence of Hawking radiation 
and discharge. A semiclassical back-reaction (T%) is considered to model Hawking radiation. We also 
included an exotic matter field ip which is a hypothetical material violating the energy condition. It is used 
to reduce the mass of the black hole substantially in a dynamic way. Discharge is simulated by supplying 
an oppositely charged scalar field to the black hole. Detailed descriptions of each process are discussed in 
the following sections. 

We use the double-null coordinate system assuming spherical symmetry 

ds 2 = —a 2 (u, v)dudv + r 2 (u, v)dQ, 2 , (3) 

where a is the lapse function, r is the radial function, and u and v indicate the in-going and out-going null 
directions, respectively [18] ■ Spherical symmetry allows a gauge choice where the gauge field is non-zero 
only in the u component: A a ~ (a(u, v), 0, 0, 0) [5]. 
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From this setup, the following components of the Einstein equation are calculated 

2 



G 



(r. uu - 2^r„) , (4) 
r V a, J 

Guv = tt^ {4rr, uv + 4r, u r,„ + a 2 ) , (5) 



2r 2 

2 / a 



— r™ — 2 — — r v ) , (6) 

r V a J 



a 

„2 



a 2 V a r a 2 / ' 

-rC 



T uu = + iea (#,«-#,«)+ eV# - ^ (8) 



47ra 2 ' 1 j 

T£ = 0,^-^, (10) 

r 2 / - - -- a v 2 \ 

Tgg = ^ ( 0>«0,t> + 0,1.0.* + * eo _ #,«) _ 2 ^,«^,D + ) ' t 11 ) 

Equations for the scalar field and the electromagnetic field are 

r<f> <uv +r, u <P,v +r, v <t>, u +iear<p tV +iear, v <t>+ -ea tV r<f> = 0, (12) 



rip,uv + r tU ip <v + r }V ip tU = 0, (13) 



i7rr 2 e 



>) = 0. (14) 



We define q = 2r 2 a tV /a 2 , which can be interpreted as the electric charge within a sphere of radius r. 

We include the semiclassical back-reaction in one-loop order from the renormalized energy-momentum 
tensor. Because of the spherical symmetry, we use the (1 + l)-dimensional result divided by 47rr 2 as an 
S-wave approximation |12)|19j 

= 4^(^-4)' < i5 > 

cs-rts = -cr^- 3 ?). < 16 > 

= 4^(^-4)' < i7 > 

where P = A^p[/127r, and AT is the number of massless scalar fields generating Hawking radiation. P 
determines the strength of the Hawking radiation. For given P, N determines the trans-Planckian scale 
where the semiclassical approximation breaks down |15j . Therefore, as we assume larger N, our approx- 
imation becomes more reliable. On the other hand, the expectation value of quantum operators has 
physical meanings on a classical background only when the quantum dispersion is relatively small. Large 
N assumption also satisfies this condition because each field independently contributes to the energy- 
momentum tensor, resulting in a small total dispersion |15j . Throughout this paper, we assume that N is 
sufficiently large. 
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For numerical integrations, we set equations following the conventions of previous researches [5][11] [ T!2 ] 



= — , h= — , / = r u , 9 = r., 
a a 



s=v47r</>, w = s >u , z = s iV , (18) 
b = \Z^rt/), o = b tU , p = b v . 



1. Einstein equations 



1 / fg a 2 a 2 q 2 1 _ _ _ _ \ 

Kv=d u = ^— + — - — -- (( wz + wz)-iea (sz-sz)) + op j : (19) 

p 

f,u = 2//i — r (ww + iea (ws ~ ws) + e 2 a 2 ss) + ro 2 (h. u — h 2 ) , (20) 

/g a 2 a 2 g 2 P 

/.» = -9-« = - — - — r + i^r - 7 d .«. ( 21 ) 

p 

= 2dg-rzz + rp 2 (cZ„-d 2 ). (22) 

r 



2. Maxwell equations 



2 

a q 
2^"' 



(23) 



3,« = (s2 - s-z) . (24) 



3. Scalar field equations 



fz gw iearz ieags ie 2 /r> \ 

= z n = ■ ■ ——a gs, (25) 

j* j> y j* /^.J*^ 

o, v = P.u = • (26) 

r r 

We solve this set of equations using the second order Runge-Kutta method [20 . h M and d )V in 
Equations (127)1) and (|2"2"j) are estimated by the finite difference method. There are three equations evolving 
the radial function r: Equations ((20"]) . (|2"Tj) . and ([22]) . Integrations using any of them should yield the same 
results. We choose Equation (f2"2"j) in most of the simulations and use the others to check the consistency 
of the code (see Appendix IX)) . 

Note that Equation (|19p is singular at r — y/P, which originates from the semiclassical approximation 
[10] . Therefore, for a reliable approximation, \[P should be sufficiently smaller than the typical length 
scale of the simulation. This constraint is satisfied throughout this paper, and, when P 7^ 0, we regard 
r = \fP as the central singularity. 

Since we are working in the double-null coordinate system, r „ = and r u = contours can be 
interpreted as the trapping and anti-trapping horizons, respectively. The electric charge q is defined above, 
and we use the Misner-Sharp mass function m(u,v) — (r/2)(l + q 2 /r 2 + 4r, M r : „/a 2 ) [21] . Therefore, M 
and Q correspond to m and q measured at the outer horizon after the black hole is formed and stabilized. 
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2.2 Initial conditions 

We simulate the gravitational collapse of a charged matter shell; we place a pulse of a complex scalar field 
along an out-going null hypersurface, and collapse it to form a black hole. We choose u = const surfaces 
as out-going null hypersurfaces and v — const surfaces as in-going null hypersurfaces. Therefore, u and v 
can be interpreted as the retarded and advanced times, respectively. We assign initial conditions for each 
function onu = and v = surfaces. 

There is a gauge freedom in choosing the initial r function. We choose r(u, 0) = urffi + and 
r(0, v) = OT,« + r^ ); therefore, g(0, v) = r,e and /(u, 0) = r,« . Because the initial condition describes a 
shell-shaped scalar field, its inside is not affected by the shell. Therefore, we can simply choose q(u, 0) = 0, 
a(u, 0) = 0, and a(u, 0) = 1. Since the mass function should vanish in this region, we choose r;« = —1/2, 



r 



( "' = 1/2 and set r<°) = r(0,0) = 10. 

We place a complex scalar field (f) along the u — surface, 



% sin2 ( cos (^TS) + ? cos (tt^ + <5)) < « < 20, 

^°' w ) = < ^^(^(cos^^+^cos^-f)) 20<«< 120, (27) 
120 < v. 



In < v < 20, a pulse with an amplitude of Ai/vAtt collapses to form a black hole. Its phase difference 5 
determines the amount of initial charge (see Equation (1241) ); field configurations with 6 — have no charge, 
and the charge increases as 5 increases until <5 = ir/2. When the discharge of the black hole is considered, 
we supply an oppositely charged pulse to the black hole, with an amplitude in 20 < v < 120. 

When the mass reduction via an exotic matter field is considered, we place tp along the u = surface, 

V>(0,0) = 0, (28) 
< v < 20, 



^ : sm (V 2=20) 20<w< 30, (29) 



v/Er V" 20 

■#= 30 < w. 

V47T 



I : i T^ v component of exotic matter field smoothly decreases from zero to —Affair) in 20 < v < 30, and 
is uniformly distributed in 30 < v. 

We place no field along the v = surface. Finally, we derive <i(0, v) — r (z~z — p 2 ) / (2g) from Equation 
(|22| assuming that the Hawking effect is negligible on the initial surfaces [12]. 

We choose the domain of computation < u < 20 and < v < 100. For all simulations, e = 0.3 and 
A\ = 0.25 are used to form an initial black hole. 
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Figure 4: Evolution of charged black holes [TS]. (a) If E + E c , discharge dominates and the black hole 
approaches E + = E c . (b) If E + <C E c , discharge is exponentially suppressed. If Q/M <C 1, Hawking 
radiation dominates and Q/M decreases, (c) If Q ^> E~ l , the black hole will approach and follow the 
extreme state, (d) If Q -C E^ 1 , the black hole will approach and follow E+ = E c . (e) When the charge 
reduces to a few quanta, Q ~ q e and Tjj ~ M _ m e on the _E + = _E C track. Then, the black hole will 
emit its final quanta of charge via Hawking radiation, (f) The final neutral black hole will lose mass via 
Hawking radiation. 
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3 Brief review of the time evolution of charged black holes 

This section summarizes the time evolution of charged black holes based on previous researches |12) |15] . 
There are two quantum effects that determine the time evolution: Hawking radiation and discharge. 
Hawking radiation is a thermal radiation of black holes, and M is reduced as the energy is emitted. It 
has Hawking temperature defined as: 



where r + = M + y/ M 2 — Q 2 is the radius of the outer horizon. The use of Misner-Sharp mass function 
guarantees that r + is equivalent to the radius of outer most trapped surface defined by r „ — 0. As the 
black hole becomes extreme, its temperature becomes zero and the radiation stops. Therefore, if one 
ignores discharge, the extreme charged black hole can be regarded as an eternal remnant. 

When the electric field near the outer horizon is strong enough, Schwinger pair-creation in this region 
reduces Q while slightly increasing M. In this regime, the pair-creation rate can be approximated as: 

q 2 E 2 e c 

I+^e-^, (31) 

where E c — irm 2 /q e is the critical pair-creating field for particles with mass m e and charge q e and 
E + = Q/r 2 ^ is the electric field near the outer horizon [35]. Note that for a given Q/M, the pair-creation 
is exponentially suppressed when the black hole is large. 

The brief paths of time evolution, considering Hawking radiation and discharge, are described in Figure 
2J which is excerpted from the previous paper |15j . For a sufficiently large black hole, pair creation is 
highly suppressed and Q/M increases via Hawking radiation. As M becomes comparable to Q, Hawking 
temperature decreases and the black hole maintains a nearly extreme state, losing its mass and charge 
simultaneously. After Q becomes sufficiently small, discharge dominates Hawking radiation and Q/M 
decreases. It keeps decreasing until the black hole emits its final quanta of charge and becomes a neutral 
black hole |23j . Finally, it will be totally evaporated, which is in the realm of quantum gravity. 



4 Formations of charged black holes with various Q/M 

Charged black holes go through states with various M and Q in their time evolution. In this section, we 
simulate gravitational collapses of matter shells with various masses and charges, considering Hawking 
radiation. The internal structures of the resulting black holes are studied, which can help in understanding 
the time evolution. 

To see the net effect of Hawking radiation, we review the static limit of charged black holes first 
[5] [11]. Figure [5] shows the result with P = 0, which is consistent with Figure [2] In a neutral black hole, 
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Figure 5: Contour diagram of radial function r and horizon structure for 6 = 0.15-7T and P = 0. This 
diagram corresponds to a dynamically formed charged black hole within the static limit. 

the singularity and the outer horizon become stuck together in the double null coordinate system (see 
Figure [6]); that is, any time-like observers beyond the outer horizon inevitably hit the central singularity. 
However, a gap arises between the singularity and the outer horizon in Figure [SJ In this region, the radial 
function monotonically decreases to some non-zero values along the out-going null direction 5 , which 
indicates the existence of the inner-Cauchy horizon at v — > oo. The curvature and the mass function 
increase exponentially in this region, which is a clear sign of mass inflation [8J . It makes the inner-Cauchy 
horizon a null curvature singularity, which preserves strong cosmic censorship. The qualitative behavior 
of charged black holes within the static limit does not vary much as we change Q/M. 

We have two free parameters that determine the amount of initial charge: e and 6. Since we are 
considering the time evolution, changing a physical constant e is irrelevant. Therefore, we change 5 to 
control the initial charge and fix the other parameters. We set an initial condition that makes a nearly 
extreme black hole when S — Tt/2, and repeat the simulation with slightly decreased 5 until 6 = and the 
black hole becomes a neutral one. We set P = 0.01 to include Hawking radiation, and A2 and A3 are set 
to zero. Among the results, we investigate some of typical ones in detail. 

4.1 5 = (Q/M = 0) 

We first discuss the 8 = (Q/M = 0) case which corresponds to an evaporating neutral black hole. The 
results are described in Figure [S] As the matter shell collapses in < v < 20, the outer horizon grows in 
a space-like direction. Then, by Hawking radiation, it bends in a time-like direction [12] |15) . Beyond the 
outer horizon, there is a space-like singularity where the curvature diverges. These results are consistent 
with the well-known characteristics of neutral black holes [231 [H] [33] . 
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6=0, P=0.01 




Figure 6: Contour diagram of radial function r and horizon structure for <5 = 0. This diagram corresponds 
to an evaporating neutral black hole. 

4.2 5 = 0.05vr (Q/M ~ 0.23) 

With 5 = 0.057T (Q/M ~ 0.23), the amount of charge is enough to generate the typical structures of 
charged black holes [12] [IS]- In Figure El r >u = and r„ = contours exist beyond the outer horizon 
that correspond to two different types of inner horizons. When Hawking radiation is not considered, an 
out-going null observer escaping from a charged black hole possibly hits the null curvature singularity 
at the inner-Cauchy horizon and then falls into the central singularity [S]. However, in an evaporating 
charged black hole, the observer reaches the r, = type inner horizon, where the radial function starts 
to increase. The r„ > region beyond this horizon causes a Cauchy horizon along the out-going null 
direction. Note that, as discussed in Section [TJ the curvature deep inside of this region becomes trans- 
Planckian because of mass inflation. Nevertheless, the horizon structure and the region around the inner 
horizon are still reliable in a background of realizable N. The finite curvature near the Cauchy horizon 
also has a possibility of violating strong cosmic censorship 15 . These regular horizons make the causal 
structure of charged black holes similar to that of nonsingular black holes |27j . 

One interesting feature is the other type of inner horizon, r u — 0. Both in a neutral black hole and a 
charged black hole within the static limit, an in-going null observer directly hits the central singularity. 
However, in an evaporating charged black hole, an r u > region exists beyond the r u = type inner 
horizon. Since this region is inside of the r \ v = type inner horizon, the radial function or, equivalently, 
the locally measured area increases for any physical observers. Therefore, as we include Hawking radiation, 
the internal structure of charged black holes becomes analogous to a throat of a wormhole. This topic is 
discussed in Section [71 
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Figure 7: Contour diagrams of radial function r and horizon structures for 8 — 0.057T. As we include 
Hawking radiation, both r tU — and r tV — types of inner horizons are observed. The resulting structure 
is similar to a throat of a wormholc. 
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6=0.15n, P=0.01 
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Figure 8: Contour diagrams of radial function r and horizon structures for S = 0.157T. As the charge is 
increased, the r u > region inside of the black hole is decreased. On the other hand, the qualitative 
behavior of the r_ v = type inner horizon is unchanged. 




Figure 9: Plot of log \M\ along u = 10 surface for black holes with 6 — 0.057T, 6\ = O.lir, and Si = 0.15-/T. 
This plot clearly shows the exponential behavior of the mass inflation effect, which becomes stronger as 
Q/M decreases. 
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Figure 10: Contour diagram of radial function r and horizon structure for S — 0.37T. Among massy 
horizons, there exists a distinct r „ = type inner horizon. The r )U > region becomes negligible. 



We investigate a black hole with more charge, 5 = 0.157T (Q/M ~ 0.61), in Figure [5] The qualitative 
behavior of the outer horizon and the r )V = type inner horizon is the same as in Figure [7] On the 
other hand, the r u = type inner horizon structure is changed from Figure [7] A boundary arises along 
an out -going null direction, where r u becomes negative again. This boundary moves downward as Q/M 
increases, and when Q/M ~ 1, the r u > region is not observed. 

The mass inflation effect also depends on Q/M of black holes. In the RN metric, mass inflation affects 
the mass function according to m(v) ~ exp(Kji;), where Ki is the surface gravity at the inner horizon 
which increases as Q /M decreases 8 . We plot the observed log \M\ in Figure |H and it clearly shows the 
exponential behavior of mass function [5] [11], Its exponent increases as Q/M decreases, which is consistent 
with the calculation in the RN metric. 



Figure [TU] describes a charged black hole with S = Q.3ir (Q/M ~ 0.91). Beyond the outer horizon, the 
contour lines of radial function become quite flat along the out-going null direction. It makes small 
fluctuations in r contour form complicated r, v = horizons, which do not seem to be very meaningful. 
However, among the messy horizons, a distinct r^ v = type inner horizon exists which approaches to the 
outer horizon. 



4.3 5 = 0.15tt (Q/M ~ 0.61) 



4.4 5 = 0.3vr (Q/M ~ 0.91) 
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Figure 11: Contour diagram of radial function r and horizon structure for 5 = 0.57T. The black hole is 
nearly extreme. Its structure is analogous to the static limit in Figure [5] 

4.5 5 = 0.5vr (Q/M ~ 0.98) 

Now we investigate a nearly extreme black hole, 8 — 0.5ir (Q/M ~ 0.98), in Figure [TT] It has very low 
Hawking temperature, which implies that the semiclassical back-reaction is negligible. One can reasonably 
expect that its structure should be analogous to the static limit in Figure [5l Interestingly, their horizon 
structures are basically the same even though their mass and charge are quite different. Neither of them 
have inner horizons nor Cauchy horizons in the scope of simulation. This observation also supports the 
physical reliability of the model. 

4.6 Summary 

The results in this section are summarized in Figure [T^l As we consider black holes with sufficient charge, 
complicated internal structures arise. Black holes with relatively small Q/M have both types of inner 
horizons, and the resulting structure is similar to a throat of a wormhole. In black holes with larger Q/M, 
the r u > region gradually disappears, and eventually, neither types of inner horizons are observed in 
nearly extreme black holes. 

It was known that Hawking radiation significantly deforms the internal structure of charged black 
holes for relatively small Q/M [HJP3]. On the other hand, Hawking radiation should be negligible for 
nearly extreme black holes. The observations in this section smoothly connect these two claims. Note 
that it cannot be applied directly to the time evolution of internal structure. This topic is discussed in 
the following section. 
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Horizon Structure Area of a Section 




Figure 12: Formation of charged black holes with various Q/M. Diagrams in the left column show the 
horizon structures and the right column show schematic diagrams of the area measured by an in-going 
null observer along the red dotted line in the corresponding left diagram. 
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5 Mass reduction via an exotic matter field 



The semiclassical back-reaction approximates Hawking radiation in a simple and proper way. Nevertheless, 
it has some drawbacks in simulating the long-term time evolution. First of all, it has an artificial singularity 
at r = \fP . Since \/P should be sufficiently smaller than the typical length scale of the simulation, 
it imposes an upper bound on the value of P. Therefore, if one wants to observe a substantial mass 
evaporation, one should carry out a long-term simulation rather than using large P for strong Hawking 
radiation. However, there is a technical difficulty that the radial contours badly fold around the asymptotic 
horizons in long-term simulations [3J . To observe the substantial mass evaporation, another approximation 
scheme is needed. 

In this section, we supply an exotic matter field to the black hole and study how the internal structures 
are affected by the substantial mass reduction. This observation is important because black holes with 
the same mass and charge can have various internal structures depending on their different histories. By 
simulating a dynamic process of mass reduction, we can get some hints on how internal structures will be 
affected during the long-term evaporation. For all simulations, we set A2 = 0, 5 — 0.157T, and P = 0. 

Figure IT51 Figure ITU and Figure [151 describe the results using A 3 = 0.001, A 3 = 0.01, and A 3 = 0.02, 
respectively. In all diagrams, time-like outer horizons and space-like inner horizons are observed, which 
is similar to Figure [S] It supports the idea of using an exotic matter field in studying the long-term 
evaporation because it affects the horizon structure in a similar way with the semiclassical back-reaction 
even though they are fundamentally different. Figure [T51 shows the Q/M of each black hole measured 
along the advanced time. As A 3 is increased, a larger amount of negative energy quickly reduces M. With 
A 3 = 0.02, Q/M reaches unity in the scope of simulation. In this case, the time-like outer horizon and 
the space-like r „ = type inner horizon meet together and disappear, which corresponds to Figure [1] (c). 
Note that Q/M cannot exceed unity effectively via Hawking radiation. 

The results above provide some implications on the long-term evaporation. Firstly, larger r u > 
region is observed in simulations with larger A3 . It is consistent with the results in Section where we 
observed larger r u > region in black holes with smaller Q/M because their high Hawking temperature 
can be matched to the large negative energy density. Secondly, the qualitative behavior of the horizon 
structure is not changed during the substantial mass reduction. This is in contrast with Section 01 where 
black holes with larger Q/M have a smaller r iU > region and no inner horizons are observed in nearly 
extreme black holes. In further simulations with various configurations of exotic matter field, we found 
that the internal structure becomes insensitive to the change of Q/M or the negative energy density with 
the lapse of time. 

To sum up, Q/M or negative energy density determines the internal structure in the stage of formation. 
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Figure 13: Contour diagrams of radial function r and horizon structures for A3 = 0.001. As the exotic 
matter field is supplied from v = 20, the outer horizon becomes time-like, and both r u = and r )U = 
types of inner horizons arise. The resulting structure is similar to the one in Figure |8l where the 
semiclassical back-reaction is used. 
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Figure 14: Contour diagrams of radial function r and horizon structures for A3 = 0.01. As a larger 
amount of negative energy is supplied, the r „ > region inside of the black hole is increased. The 
relation between the negative energy and the wormhole-like structure is discussed in Section [7] The outer 
horizon approaches to the inner horizon more quickly. 
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6=0.15n, P=0,Aj=0.()2 




Figure 15: Contour diagrams of radial function r and horizon structures for A3 = 0.02. As Q/M reaches 
unity, the outer horizon and the r v = type inner horizon meet together and disappear, which is consistent 
with the RN metric in Figure [T] (c). On the other hand, the r „ = type inner horizon maintains its 
structure. 




Figure 16: Q/M measured at the outer horizon along the advanced time v. 
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Figure 17: Long-term evaporation of charged black holes. Diagrams in the lower row show schematic 
diagrams of the area measured by an in-going null observer along the red dotted lines in the upper 
diagram. 

However, it becomes insensitive to them with the lapse of time, and the structure is fixed. The reason 
seems to come from mass inflation. It induces large curvature to the internal structure, which makes it 
insensitive to the late-time perturbations. Eventually, the curvature will become trans-Planckian. In this 
regime, the inner horizon can be regarded as a physical singularity with a finite radius, which is analogous 
to some dilaton black holes models [15] [28] . 

We summarize the results in this section in Figure 1171 

6 Neutralization of charged black holes 

As a charged black hole gets some opposite charge, it discharges and becomes a neutral black hole. In this 
section, we investigate the neutralization of charged black holes for both non-evaporating and evaporating 
cases. For all simulations, we set = 0.055, 8 = 0.157T, and P = 0. 

We first investigate the neutralization of non-evaporating charged black holes, described in Figure IT51 
Compare this result with Figure[5] which describes the same initial black hole, but not the neutralized one. 
As the black hole discharges, the radial function beyond the outer horizon decreases, and the singularity 
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6=0.15n, P=0, As=0.055,As=0 




Figure 18: Contour diagram of radial function r and horizon structure for A3 = 0. The charged black 
hole in Figure [5] discharges and becomes a neutral black hole. 

moves downward. As a result, the gap between the singularity and the outer horizon gradually disappears, 
and the structure becomes a neutral one. An out-going null observer who would reach the inner-Cauchy 
horizon in Figure [5] will hit the central singularity if the black hole is neutralized. 

In the previous paper, the neutralization of evaporating charged black holes is studied using the S- 
wave approximation of semiclassical back-reaction |15] . The key feature of this process is the inner horizon 
evolving to a singularity. However, the S-wave approximation causes an artificial singularity that disturbs 
a direct observation of the formation of the singularity. To avoid this problem, we adopt an exotic matter 
field to produce an evaporating internal structure and re-test the neutralization process. We add an exotic 
matter field with A3 = 0.01 to the discharging black hole in Figure [TH] The result is displayed in Figure 
1191 Similar to the non-evaporating case, the radial function beyond the outer horizon decreases as the 
black hole is neutralized. However, the local minimum of radial function for a given v is the r jV = 
type inner horizon in evaporating charged black holes. Therefore, the decreasing radial function becomes 
a singularity at the inner horizon, and the region beyond the inner horizon is causally disconnected, 
leaving the structure of neutral black holes. This result is the same as the analysis using the semiclassical 
back-reaction, which confirms the previous results |15j . 

We summarize the neutralization of evaporating charged black holes in Figure 

7 Charged black hole as a wormhole— like structure 

The RN metric has wormholes called RN wormholes [T] [5] . However, they are not observable in dynamically 
formed charged black holes within the static limit because of the mass inflation effect. The wormhole-like 
structure appears again as we consider Hawking radiation. Evaporating charged black holes with relatively 
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Figure 19: Contour diagrams of radial function r and horizon structures for A3 = 0.01. An exotic matter 
field is added to the discharging black hole in Figure [THI to consider Hawking radiation. 
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Figure 20: Neutralization of evaporating charged black holes. Diagrams in the lower row show schematic 
diagrams of the area measured by an in-going null observer along the red dotted lines in the upper 
diagram. 
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Figure 21: The neutralization of charged black holes can be identified with the separation between the 
inside and the outside of the wormhole-like structure, (a) If the inside region is regular and does not have 
a time-like singularity, this process can cause issues related to the information loss problem, (b) In other 
cases, semiclassical analysis cannot be applied in this region. 
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small Q/M have both r „ = and r v — types of inner horizons, and any physical observers crossing 
both of them will observe the increase of locally measured area. The inner horizon corresponds to the 
throat of the wormhole-like structure, as described in Figure [T2l 
In Equation G uu component of the Einstein tensor gives 

- - (r,uu - 2— r„] = 8nT uu , (32) 

where T uu = T^ u + (T^ u ). To obtain a wormhole-like structure, r^ uu > around the r M = type horizon 
is necessary. It requires T uu < in this region, violation of the null energy condition. Since we assumed 
Hawking radiation, such violation is naturally obtained. In contrast to many models of wormholes that 
involve exotic matters or exotic initial conditions, the wormhole-like structure in evaporating charged 
black holes can be an interesting and natural model to obtain such a structure in a dynamic process 
29J [30 31 [22] [33]. However, it is not clear whether a charged black hole is indeed a wormhole or not. 
While a wormhole connects two asymptotically fiat regions, the existence of such a region inside of a 
charged black hole cannot be determined by semiclassical analysis because of the Cauchy horizon. Even 
if a charged black hole is a wormhole, it is one-way traversable and its deep inside has large curvature 
which can be trans-Planckian depending on N. Nevertheless, the region around the throat is still reliable 
in a background of realizable N. 

As a charged black hole is neutralized, its inner horizon evolves to a singularity. During this process, the 
inside and the outside of the wormhole-like structure are separated and the throat becomes a singularity. 
It induces a Cauchy horizon, beyond which cannot be determined in terms of the semiclassical physics. 
Nevertheless, suppose that the inside region is regular and that there is no time-like singularity, as 
described in Figured]] (a). In this case, this process can cause issues related to the information loss problem 
because the disconnected region can have an asymptotically flat region 15 [34 [35 36 [37]. However, in 
other cases such as Figure [21] (b), the semiclassical analysis does not provide any meaningful remarks. 

8 Conclusion 

We studied the internal structure of charged black holes for important stages in its time evolution using 
numerical simulations. We adopted a spherically symmetric charged black hole model that includes the 
renormalized energy-momentum tensor to consider Hawking radiation. Solving the equations for various 
initial conditions, we obtained general behavior of charged black holes. 

When charged black holes are formed by collapses of charged matter shells, various internal structures 
are observed depending on their mass and charge. For black holes with relatively small charge-mass ratio, 
both trapping and anti-trapping inner horizons are observed. The internal structure of such black holes 
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is similar to a wormhole. As we consider black holes with more charge, the inner horizons gradually move 
out of sight, and they are not observed in nearly extreme black holes. Since Hawking temperature is close 
to zero in nearly extreme black holes, it is natural that their structure becomes analogous to the static 
limit. 

As time passes, the mass and charge of the black hole will evolve via Hawking radiation and discharge. 
Although the mass and charge determined the internal structure in the stage of formation, we found that 
the internal structure becomes insensitive to them with the lapse of time. The reason seems to come from 
mass inflation which induces large curvature in the internal structure. The internal structure determined 
from the formation seems to be maintained during the evaporation. 

We also investigated the neutralization of charged black holes for both non-evaporating and evapo- 
rating cases. As charged black holes obtain some opposite charge, the radial function beyond the outer 
horizon decreases in both cases, causing the singularity to approach to the outer horizon rapidly in non- 
evaporating charged black holes. In evaporating charged black holes, obtaining of opposite charge causes 
the inner horizon to evolve to a singularity. 

Finally, we focused on the wormhole-like structure inside of charged black holes which is a worthwhile 
object of examination because, without any exotic requirements, Hawking radiation naturally induces 
such a structure. From this point of view, the neutralization of charged black holes can be interpreted 
as the separation between the inside and the outside of the wormhole-like structure. It may cause some 
issues related to the information loss problem; however, further study is required. 
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A Convergence and consistency checks 

In this appendix, we check the convergence and consistency of the numerical calculations. The initial con- 
dition described in Figure [5] was tested in this section. We also successfully tested other initial conditions, 
which include non-zero A2 and ^.3. 
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Figure 22: |rm — ^(2)1/^(2) & n d M r (2) ~ r (4) I / r (4) along a few constant u lines, where r( n ) is the radial 
function calculated innxn times finer grid than rm. It converges to the second order with errors < 0.1%. 




Figure 23: Differences in radial function r between two integration schemes along a few constant u lines. 
r(„) is calculated by integrating g using Equation (|2"2")l . and r( u ) is calculated by integrating / using 
Equation (gD]). The error is < 1%. 
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To check the convergence, we compared the results from different step sizes: 1 times, 2 times, and 4 
times finer for each axis. Figure l2"2l shows that the differences in radial function r decrease by 1/4 as the 
step size becomes half. The differences converge to the second order with errors < 0.1%. 

To check the consistency, we compared the results from two different sets of equations. In most of the 
simulations, we used Equations (|2ip and (1221) to evolve / and g, respectively, and d, g, and z to evolve a, 
r, and s, respectively. We compare this result with the other one using Equations (l20l) and (|21[) to evolve 
/ and <?, respectively, and h, /, and w to evolve a, r, and s, respectively. The radial function r integrated 
from each scheme corresponds to rr v \ and rr u ) in Figure [2"31 respectively. The difference is < 1%. 

Near the horizons and the singularity, the errors increase rapidly, and it comes from the rapidly varying 
radial function in these regions. We adopted the adaptive grid refinement method (AGM) to handle this 
problem [5]. In this paper, we adaptively set the step size Au in u direction so that r tU Au/r < 1%. 
When the radial function varies rapidly, AGM automatically reduces the step size so that more detailed 
calculations can be done in that region. 
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